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Non-Fickian Diffusion and the Accumulation of Methane 
Bubbles in Deep- Water Sediments 
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In the absence of fractures, methane bubbles in deep- 
water sediments are immovably trapped within a porous 
matrix by surface tension. The dominant mechanism of 
transfer of gas mass therefore becomes the diffusion of gas 
molecules through porewater. The accurate description of 
this process requires non-Fickian diffusion to be accounted 
for, including both thermodiffusion and gravitational action. 
We evaluate the diffusive flux of aqueous methane consider- 
ing non-Fickian diffusion and predict the existence of bubble 
mass accumulation zones within deep-water sediments. The 
occurrence of these may be highly significant in the assess- 
ment of methane hydrate reservoirs or inventories as they 
could occur independently of the hydrate stability zone, yet 
may mimic the bottom-simulating-reflector which is com- 
monly used to identify the base of the zone. 



1. Introduction 

The occurrence of methane bubbles within porous water- 
saturated sediments is widespread around the ocean mar- 
gins. The gas within them plays an important role in both 
submarine hazards, such as submarine landslides [Kayen and 
Lee, 1991] as well as the formation of resources, such as 
methane-hydrate deposits [Davie and Buffett, 2001, 2003; 
Archer, 2007]. The stability of the bubbles has a signifi- 
cant control on the methane flux from the sediments into 
the ocean-atmosphere system. 

In porous sediments the bubbles are trapped within the 
matrix pores. Large moving bubbles are unstable, as they 
split into smaller bubbles during migration [Lyubimov et al, 
2009; Barry et al., 2010], and smaller bubbles are trapped 
by pore-throats or by surface tension forces. The mini- 
mum pore throat diameter I required to trap a small bub- 
ble, when there is no strong pumping of the fluid through 
the porous matrix, can be calculated. The surface ten- 
sion forces al (a is the surface tension) should overwhelm 
the buoyancy force pgl 3 (p is the density, g is the grav- 
ity), i.e., / < \Ja/ pg. For gas- water systems, one finds, 
I < [0.073 N m" 1 /10 3 kg m,- 3 9.8 ms~ 2 ] 1/2 w 2.7mm. Mak- 



1 Department of Mathematics, University of Leicester, 
Leicester LEI 7RH, UK. 

2 Department of Geology, University of Leicester, Leicester 
LEI 7RH, UK. 

3 Institute of Continuous Media Mechanics, UB RAS, 
Perm 614013, Russia. 

4 British Geological Survey, Keyworth, Nottingham NG12 
5GG, UK. 

5 School of Earth and Environment, University of Leeds, 
Leeds LS2 9JT, UK. 



Copyright 2011 by the American Geophysical Union. 
0094-8276/ll/$5.00 



ing allowance for the inhomogeneity of pores and the ge- 
ometry of contacts, one should decrease this estimate to 
I ~ 1 mm, which still suggests trapping even for sands. For 
a soft mud, mechanics of bubbles and porous matrix can be 
different [Algar et al., 2011] and is not considered here. 

As bubbles are trapped, and in the absence of significant 
groundwater movement transporting dissolved methane, the 
dominant mechanism of methane mass transfer in deep- 
water sediment is by diffusion of methane through porewa- 
ter, controlled by (i) the methane saturation of the aque- 
ous solution throughout the sediment volume, and (ii) non- 
Fickian diffusion laws. The latter pertain to two processes: 
firstly, the geothermal gradient causes thermodiffusion (the 
Soret effect [Soret, 1879]), where the temperature gradi- 
ent induces solute flux (as recognized in several fields, e.g. 
[Richter, 1972]); and secondly, the impact of gravity on dis- 
solved molecules. These non-Fickian contributions to the 
diffusive flux mean that the solute flux cannot solely be de- 
termined by the gradient of the solute concentration. 

Currently non-Fickian diffusion is rarely considered in the 
modelling of deep-water sediment systems (Davie and Buf- 
fett [2001, 2003]; Haacke et al. [2008]; Garg et al. [2008], 
and others involved in gas hydrate modelling only address 
Fickian diffusion). However, we suggest that non-Fickian 
processes are important; whilst appearing counter-intuitive, 
they may cause methane to migrate against the direction 
of the steepest decrease of concentration under certain con- 
ditions. Unfortunately, the experimental value of the ther- 
modiffusion coefficient for the aqueous methane solution is 
unknown and can only be roughly assessed theoretically. 
Therefore, we treat it as a free parameter in our investi- 
gation. This uncertainty also justifies the fact that prior 
numerical modelling work did not include non-Fickian pro- 
cesses. 

In researching the horizontal through-flow of water and 
vertical aqueous oxygen transport in porous bubble-bearing 
sediments, Donaldson et al. [1997] considered the hydro- 
dynamic dispersion (or "turbulent diffusion" [Bernard and 
Wtlhelm, 1950]), caused by water transport through irregu- 
lar pore channels, as a transport mechanism and reasonably 
neglected the molecular diffusion. Indeed, the turbulent dif- 
fusion plays a significant role in vertical dissolved gas trans- 
port near the earth's surface. However, it becomes insignif- 
icant away from the sediment-water interface in deep-water 
sediments, where vertical and horizontal displacements of 
water are comparable. Here it is significantly smaller than 
molecular diffusion — even in sandy sediments. 

In the present study we consider the diffusive migration 
of methane in sea-floor sediments where water is saturated 
with methane, and some methane is gaseous (forming bub- 
bles; see sketches in Fig. 1). To evaluate the solubility and 
diffusive flux of methane we employ the physical model de- 
veloped by Goldobin and Brilliantov [2011] and summarized 
in the Appendix. We demonstrate that even in the presence 
of the advective transport by the upward fluid flux, the non- 
Fickain diffusion is still important. 
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Figure 1. Four possible profiles of diffusive flux driven 
by the solubility profile and non-Fickian effects. The blue 
solid lines depict the flux J vs the depth below the water- 
sediment interface, z. The red dashed lines correspond 
to the free-gas production Xb, resulted from the diffusive 
flux J. In the diagrams the directions of fluxes and zones 
of methane accumulation from aqueous solution into free- 
gas bubbles (that is zones of positive Xb) are shown. 



2. Diffusion in non-isothermal 
solutions bearing bubbles 



aqueous 



Under non-isothermal conditions the diffusive flux of so- 
lute molar fraction X in the solvent is governed by the law 
(cf [Bird et al., 2007; Goldobin and Brilliantov, 2011]) 



Jdiff 



-<j)DX 



X7X VT 
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Mg 
RT 
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Here D is the solute effective molecular diffusion coefficient, 
cj> is the porosity of the solid matrix. The first term de- 
scribes the "ordinary" Fickian diffusion, JFi c k = -~4>DV X. 
The second term represents the thermodiffusion effect ap- 
pearing in non-isothermal systems, where temperature in- 
homogeneity causes a solute flux. The strength of the ther- 
modiffusion effect is characterized by the thermodiffusion 
constant a (the conventional Soret or separation coefficient 
St = oe/T). The third term describes the action of gravity 
on solute molecules; R = 8.31 J / (mol K) is the universal gas 
constant, M = M B - AiA/ host , M g and M hoBt are the molar 
masses of the solute and solvent, respectively, and iVi is the 
number of solvent molecules in the volume occupied by one 
solute molecule in the solution. The value of Ni can be pre- 
cisely derived for CH4-H2O systems from the dependence of 
the solution density on its concentration [Hnedkovsky et al, 
1996]; one obtains iVi = 2.23 and M = -24.3 g/mol. 

When the liquid is saturated with gas bubbles, the con- 
centration of solute in solvent equals the solubility, X — 
X^°\ throughout the liquid volume because the bubbles 
are in local thermodynamic equilibrium with the solution. 
Thus, the solute flux depends merely on the temperature 
and pressure fields, T(z) and P(z), and the solution concen- 
tration is not a free variable, X(z) = X (0) (T(z), P{z)). For 
the calculation of the solubility see the Appendix. 

In marine sediments, the processes of sediment com- 
paction and related fluid advection can be important [Davie 



and Buffett, 2001, 
flux reads 



2003]. With these processes, the solute 



J = Jdiff + UfX, 



(2) 



where Uf is the fluid filtration velocity, which is related to 
the interstitial fluid velocity vt = ut/4>(z). 

Generally, the flux (2) possesses non-zero divergence 
V • J, which requires sources and sinks of the solute mass — 
provided by bubbles. The mass production ( — V ■ J) is con- 
sumed by the bubble volume; 



dX b 
dt 



+ v s • VX b 



= -V- J, 



(3) 



where Xb is the molar fraction of bubbles in the bubble- 
bearing fluid (comprising bubbles and liquid); v s is the ve- 
locity of sediments movement due to the compaction. By 
virtue of the mass conservation law for sediments, v s (z) = 
([1 - 0(O)]/[1 - (t>(z)])v s (0) [Davie and Buffett, 2001]. 

3. Evolution of methane bubbles in the 
seabed 

On the field scale, deep-water sediments are typically 
much more uniform horizontally rather than vertically. Con- 
sequently, we consider a system that is uniform horizon- 
tally. The depth below the seafloor is measured by the z- 
coordinate (Fig. 1). The system is essentially characterized 
by the hydrostatic pressure P and geothermal temperature 
gradient G; 



P(z) = P +pi icl g(z + H), T(z)=T s{ + Gz, (4) 

where Po is the atmospheric pressure, H is the height of 
sea-water (water depth) above the bubble-bearing porous 
sediments, and T B f is the temperature at the sediment- water 
interface. 

For given hydrostatic pressure and geothermal gradient, 
the solute flux (1) reads: 



Jdiff = <j>D 



where 



dz 



z + T st /G 



(5) 



Mg 
RG' 



(6) 



Eq. (5) with /3 = corresponds to the Fickian diffusion law, 
^Fick = —<f)DX7X, and /3 characterizes the strength of non- 
Fickian part of the solute flux. As the value of a is unknown, 
/3 is treated as a free parameter in our study. 

3.1. Diffusive transport 

First, we focus on purely diffusive transport without 
advection, which is reasonable for the following reasons: 
Firstly, for some geological systems the molecular diffusion is 
expected to be the dominant transport mechanism [Haacke 
et al., 2008]. Secondly, advection is a well studied process; it 
is not influenced by the non-Fickian effects addressed here. 
In the saturated solution (the bubbly zone) advective and 
diffusive transports are merely summed. 
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Figure 2. Diagram of diffusive methane transport regimes on the fi-H plane (/3 quantifies the non-Fickian drift strength), 
in the seabed with methane bubbles, for the water-sediment interface temperature T s f = 277 K. The geothermal gradient 
is G — 40 K /km (a) and G — 60 K /km (b). The profiles of solute flux J for 4 possible regimes (A, B, C, D) are plotted in 
Fig. 1. 



We perform the following calculation: Employing Eq. (7) 
from the Appendix, we find the water-free-gas solubility 
profile for the temperature and pressure profiles given by 
Eq. (4) (such kind solubility profiles can be found in [Haacke 
et ai, 2008]). For the calculated solubility profile X (0) (z) 
and parameter /3, Eq. (5) yields the profile of the diffusive 
flux J(z). The variety of the solubility profiles for possible 
values of the temperature of the sediment-water interface 
T s f and the depth of the water body above sediments H in 
combination with various strengths of non-Fickian drift ft 
result in 4 classes of the flux profiles plotted in Fig. I. Ac- 
cording to Eq. (3), these flux profiles cause the production 
of free-gas mass Xb = —(dj/dz) for Xb > or its depletion 
for X h < 0. 

We now discuss the diagrams of these regimes on the /3- 
H plane. Fig. 2a shows the diagram of diffusive regimes in 
the system with geothermal gradient G = 40 K/km. For 
shallow water bodies (small H), regime (A) occurs, with an 
upward diffusive flux near to the sediment-water interface, 
but a downward flux deeper within the sediments. As water 
depth increases, the flux inversion point shifts towards the 
sediment-water interface and finally disappears, so that the 
zone of downward diffusive methane flux covers the whole 
sediment volume, as shown in Fig. 1, regime (B). In regimes 
(A) and (B) all methane leaves the upper 2 km layer of sed- 
iments either upwards, as identified locally in (A), or down- 
wards, as in (B), and locally in (A). Intriguingly, besides 
regimes (A) and (B), there are regions where methane ac- 
cumulation zones occur, caused by both upward diffusion 
(from below) and downward diffusion (from above). These 
zones may occur within the sediment column (C) and can 
touch the sediment- water interface (D). For G = 40 K/km 
this methane accumulation zone exists for non-Fickian drift 
strength j3 in range [0.2, 1.1] in water bodies with depth be- 
tween 0.5 km and 1.5 km, and /3 < —0.68 in water bodies 
with depths greater than 3.2 km. 

The regime diagram in Fig. 2a is plotted for a geother- 
mal gradient G = 40 K/km, which is typical for the Blake 
Ridge [Paull et ai, 2000; Davie and Buffett, 2001]. As G 



increases, regions (C) and (D), and the boundary between 
regions (A) and (B), shift to smaller depths and negative 
values of non-Fickian drift strength /3. Fig. 2b presents the 
regime diagram for G — 60 K/km, which is typical for the 
Cascadia margin [Davie and Buffett, 2003]. At sites with a 
low geothermal gradient, G = 20 K/km, regime (A) prevails 
for water depths up to 5 km (not illustrated). 

Note, that the diagrams in Fig. 2 are constructed 
analysing the direction of flux J^is within the sediment col- 
umn. The sign of mass production Xb is also affected by the 
nonuniformity of <f>(z) (see Eqs. (3) and (5)). Importantly, 
the positive mass production zones, Xb > 0, in regimes (C) 
and (D) appear for any dependence <f>{z). In the following 
subsection we calculate Xb, taking into account the pres- 
ence of an upward fluid flow, as well as the nonuniformity 
oi(j>{z). 

3.2. Impact of advection 

Advective transport is more system-specific than diffu- 
sive flux; we need to additionally specify the dependence of 
the porosity <f) on the vertical coordinate z, the sedimenta- 
tion rate v B (0), and the fluid filtration velocity Uf. For the 
west Svalbard continental slope Haacke et al. [2008] adopted 
the dependence cj>(z) = 0(0) exp( — z/L) with (f)(0) — 0.555 
and L = 1053m. Sedimentation rate u a (0) = 50 cm/kyr, 
geothermal gradient G = 86.5 K /km, molecular diffusion 
coefficient D — 2 ■ 10~ 9 m 2 /s, and the upward fluid velocity 
assessed as ff(0) = 0.1 mm/yr have been used. Employ- 
ing these parameters we calculated from Eq. (3) the free-gas 
mass production from the saturated aqueous solution. For a 
better understanding of the impact of advection the calcula- 
tions have been performed not only for V{ (0) = 0.1 mm/yr, 
but also for smaller values of the upward fluid flux. 

In Fig. 3 we plot the zones where methane is accumu- 
lated from the saturated solution into free-gas bubbles and 
the zones where free gas is dissolved. The lines demarcate 
the accumulation and dissolution zones for specified values 
of non-Fickian drift strength (3. The former are observed 
for higher depths of water bodies next to the gas-hydrate 
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Figure 3. Boundaries between the zones of the growth of free-gas bubbles from the aqueous solution and the zones of 
dissolution of bubbles for different strength of non-Fickian drift f3. The upward fluid flux Uf(0) = Uf/(j)(0) = 0.02mm/yr 
(a), Vf(0) = 0.05 mm/yr (b), and Vf(0) = O.lmm/yr (c). 



stability zone (GHSZ), and the latter occur under shallow 
water bodies. We wish to stress that the bubble growth (or 
dissolution) in the bubbly zone is a local process, which is 
not affected by the neighboring GHSZ (if any). This allows 
one to exclude the GHSZ from the explicit consideration. 
One can see that for a stronger advective flux Vf, the ac- 
cumulation zone shrinks and its existence requires bigger 
positive /3. As one can see in Fig. 3, the non-Fickian correc- 
tion to molecular diffusion remains rather important in the 
presence of advection. 

3.3. Importance and uncertainties 

The behavior described is considerably influenced by the 
non-Fickian drift of methane. We suggest that the Fickian 
diffusion law, which has been adopted in [Davie and Buffett, 
2001] and its successors and corresponds to /3 = 0, should 
be modified. 

The unresolved issue here is the specific value of /3 for 
methane. The authors are not aware of experimental data on 
thermodiffusion of methane in water, though there are a lot 
of experimental studies on the thermodiffusion of methane in 
mixtures of hydrocarbons (e.g. [Wittko and Kohler, 2005]). 
Theoretical studies (e.g. [Semenov, 2010]) provide formulae 
for calculation of the thermodiffusion constant from inter- 
molecular potentials which are poorly established for wa- 
ter because of hydrogen bonds. We can only calculate 
Mg/RG = -0.725 (for G = 40 if/fern), the isotope (or ki- 
netic) part of the thermodiffusion [Wittko and Kohler, 2005; 
Semenov, 2010] a iso t = (3/4) ln(Af g /iViM host ) = -0.670, 
and infer that the inter-molecule potential part a po tcnt is 
positive (this qualitative conclusion is suggested by formu- 
lae in [Semenov, 2010]). Thus, /3 CH4 < -0.055 for aque- 
ous solutions with geothermal gradient G = 40 K/km, and 
/3 CH4 < 0.187 for G = 60 K/km. 

4. Implications for marine hydrates 

The processes described above have major significance in 
relation to the occurrence of gas hydrate deposits. They sug- 
gest that in relatively shallow waters (area (A) in Fig. 2) the 
diffusive methane flux is directed upwards, into the GHSZ, 
while in deeper waters (area (B) in Fig. 2) methane dif- 
fuses downwards, leaving the GHSZ. In the former case the 
hydrate deposit capacity is enhanced, while in the latter, 
the deposit suffers diffusive depletion. The threshold sea 
depth, where the transition between these two regimes oc- 
curs, strongly depends on the non-Fickian drift strength ft 
(which is a function of the thermodiffusion constant a and 



the geothermal gradient G). Nevertheless, theoretical as- 
sessment of the value of a suggests enhancements of hydrate 
deposits for the Cascadia margin and the Blake Ridge. 

Another implication relates to the detection of marine 
methane hydrates. A very small amount (c. 1 — 2%) of gas 
bubbles is enough to change the sound speed in the sedi- 
ment body and cause reflection of seismic waves; the reflect- 
ing layer is commonly referred to as the 'bottom simulating 
reflector' (BSR) (e.g. [MacKay et al, 1994; Hovland et al, 
1999; Paull et al., 2000]). It is typically assumed that the 
BSR marks the base of the GHSZ. Indeed, for an isothermal 
Fickian diffusion law, there are no reasons for the forma- 
tion of a gaseous methane horizon (that may cut-across pri- 
mary sedimentary features, such as bedding) other than the 
base of the GHSZ above which gas disappears having been 
transformed into hydrate. However, our theoretical findings 
suggest that commonly observed geothermal gradients are 
enough to explain the creation of gaseous methane horizons 
by thermodiffusion; formation of these horizons does not re- 
quire the presence of the GHSZ or hydrates therein. This 
may be a reason why some BSRs seem to have no hydrate 
associated with them (e.g. [Paull et al., 2000]). Rigorous 
quantitative prediction of appearance of BSRs without hy- 
drates requires consideration of hydrate stability and will be 
presented in a separate publication. With the present work 
we conclude that within regions (C) and (D) in Fig. 2 con- 
ditions are favorable for the creation of potential BSRs that 
are generated by thermodiffusion processes alone, and are 
not linked to the boundary of a GHSZ at that depth. The 
Cascadia margin with H — 1.3 km and G = 60 K/km may 
well be close to regime (D) (Fig. 2b). 

Our results suggest that some hydrate deposits may have 
no BSR at their base (as, for instance, reported by Ecker 
et al. [1998]; Paull et al. [2000]) because diffusion can lead 
to depletion of methane from the bubbly area as it does in 
regimes (A) and (B), which are most common (Fig. 2). In 
this case, only the balance between the sedimentation rate 
(which provides a flux of methane mass from the GHSZ to 
the underlying free-gas zone) and the diffusive depletion rate 
determines whether the BSR appears. In areas with higher 
sedimentation rates it may be that BSRs will be more com- 
mon. This relationship should be addressed in future mod- 
elling, involving consideration of methane-hydrate stability 
and sedimentation processes (as in [Davie and Buffett, 2001], 
but taking account of the non-Fickian effects in diffusion) . 

5. Conclusion 

We have theoretically explored the process of diffusive 
migration of aqueous methane in the presence of bubbles, 
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when they are immovably trapped by a porous matrix — 
as occurs commonly in seafloor sediments, swamps, or ter- 
restrial aquifers. The effect of temperature inhomogeneity 
across the system (geothermal gradient) and gravitational 
force have been accounted for. 

Non-Fickian corrections — thermodiffusion and gravita- 
tional buoyancy — appear to play an important role in the 
migration of methane in sediments in deep-water settings. 
They can assist the formation of a methane accumulation 
zone in the upper part of the sediment column under water 
bodies with depths in range 0.5 — 1.5 km and the very deep 
ones (Fig. 2). For the latter the minimal depth required de- 
pends on the geothermal gradient G: it is around 3.2 km 
for G = 40 K/km (such G is typical for the Blake Ridge) 
and 2.2 km for G = 60 K/km (such a G is typical for the 
Cascadia margin). 

Under moderately deep water bodies, gas hydrate de- 
posits should be enhanced (presumably, as for the Casca- 
dia margin and the Blake Ridge), and conversely, under 
deep water bodies gas hydrate deposits should be subject 
to a stronger diffusive depletion than with the purely Fick- 
ian diffusion. The depth of transition between these two 
regimes increases with decrease of the geothermal gradient 
G: this is the reason why the relatively deep Blake Ridge 
with G « 40 K/km and the shallower Cascadia margin with 
large G « 60 K/km appear to have similar conditions for 
the formation of gaseous methane accumulation zones. 

It is a consequence of our results, that the effects of non- 
Fickian diffusion remain rather important even in the pres- 
ence of advective transport. For instance, Fig. 3 illustrates 
that, depending on f3 and for conditions of the west Svalbard 
continental slope [Haacke et at, 2008], non-Fickian diffusion 
can either extend or shrink the zone of methane accumula- 
tion from solution into bubbles. 

Remarkably, the early theoretical study [Xu and Rup- 
pel, 1999] on the formation of natural methane hydrate de- 
posits had inferred that these systems should be advection 
dominated; this was needed to explain the enhancement of 
hydrate deposit capacity. Later on, Haacke et al. [2008] 
reported a case study from the west Svalbard continental 
slope, where the formation of a persistent free-gas zone be- 
low the gas hydrate stability zone required the system to 
be diffusion dominated. These studies have adopted the 
Fickian diffusion law. In our study we see that the actual, 
non-Fickian diffusion, even for negligible advection, can lead 
both to the enhancement of hydrate deposits and the for- 
mation of a persistent free-gas zone. 

Unfortunately, we cannot determine precise values for the 
thermodiffusion constant of aqueous solutions of methane 
from the literature, and can only rely on theoretical predic- 
tions (e.g. [Semenov, 2010]), to estimate their values, as we 
do here. Our findings highlight the necessity of experimental 
determinations of the thermodiffusion constant for aqueous 
methane solutions. 
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Appendix: Methane solubility 

For theoretical evaluation of the methane solubility at 
equilibrium with the vapor phase, we employ the scaled par- 
ticle theory [Pierotti, 1976] and the Van der Waals' real gas 
model (see [Goldobin and Brilliantov, 2011] for mathemati- 
cal details). These yield 



y(0) _ (1 - Y)v liq n 
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Here G c is the enthalpy of formation of the cavities for one 
mole of solute molecules, 
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the effective diameter of water molecules <7ii q = 2.77 A, the 
one for methane a g = 3.27 A; y = (7r/6)nii q cr 1 ; q (at standard 
conditions, yn 2 o = 0.371, the corresponding molar volume 
of liquid water t>i; q = 17.95 cm 3 /mol [Pierotti, 1976]); Ava- 
gadro's number N A = 6.02- 10 23 . The molar density n(T, P) 
of gaseous methane can be determined from the Van der 
Waals equation of state 

P = nRT/(l - nb) - an 

with a = 0.2283 m 6 Pa /mol 2 and b = 4.278 ■ W~ 5 m 3 /mol, 
which has unique analytical solution with respect to n for 
T > 190 K (the expression for this solution is lengthy though 
can be provided by any package for analytical computa- 
tions, such as Maple). The molar enthalpy of interaction 
of methane molecule with the surrounding water molecules, 
Gi, is nearly temperature independent, Gi/R ~ — 1138.K\ 
The molar fraction of water in the vapor phase 
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°(--—\ 
\T To) 



(8) 



where the difference between the molar specific heat capac- 



ity of liquid water and steam Ac p 2 / R = —5.00, the molar 
enthalpy of vaporization Ah% 2 ° /R = 4890 K at T = 373 K 
and Po = 1 atm. Thus, one has all ingredients for calcula- 
tion of the solubility A (0) (T, P) with Eq. (7). 



Notation 

a, b Van der Waals' constants, Pam 6 /mol 2 and 
m 3 /mol, respectively. 

D molecular diffusion coefficient, m 2 /s. 
G geothermal gradient, K/m. 
g gravity, m/s 2 . 
G c enthalpy of formation of the cavities in the so- 
lution for one mole of the solute molecules, 
J/mol. 

d enthalpy of interaction of the solute molecule 
with the solvent molecules, per one mole of the 
former, J/mol. 

H depth of the water body above sediments, m. 

J diffusive flux of the CH4 molar fraction, m/s. 

M effective molar mass of methane, M — Mch 4 — 
A^iMh 2 o, where N\ is the number of the water 
molecules replaced by the methane molecule in 
the solution; M = -0.0243 kg/mol. 
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n the molar density of gaseous methane, mol/m 3 . 

P, Pq pressure and the atmospheric pressure, respec- 
tively, Pa. 

R universal gas constant, 8.31 J/(mol K). 
T, T s f temperature and the sediment-water interface 
temperature, respectively, K. 
Uf filtration velocity of fluid, m/s. 
Vf interstitial fluid velocity, Uf = (jwf, m/s. 
viiq molar volume of liquid water, m 3 /mol. 

X, X (0) molar fraction of CH 4 in aqueous solution and 
the solubility, respectively, dimensionless. 
At, molar fraction of gaseous methane in bubbly 
fluid, dimensionless. 
Y molar fraction of water in the methane bubble, 

dimensionless. 
z depth in sediments, below the sediment-water 
interface, m. 

q thermodiffusion constant, dimensionless; Soret 
coefficient St = a/T. 

(3 — —a + Mg/RG characterizes the non-Fickian 
drift strength, dimensionless. 
<p(z) porosity of sediments, dimensionless. 
piiq density of liquid water, kg/m 3 . 
BSR bottom simulating reflector. 
GHSZ gas hydrate stability zone. 
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